Estimating the costs and perceived benefits of oral pre-exposure prophylaxis (PrEP) delivery in ten counties of Kenya: a costing and a contingent valuation study.

Background: Kenya included oral PrEP in the national guidelines as part of combination HIV prevention, and subsequently began providing PrEP to individuals who are at elevated risk of HIV infection in 2017. However, as scale-up continued, there was a recognized gap in knowledge on the cost of delivering oral PrEP. This gap limited the ability of the Government of Kenya to budget for its PrEP scale-up and to evaluate PrEP relative to other HIV prevention strategies. The following study calculated the actual costs of oral PrEP scale-up as it was being delivered in ten counties in Kenya. This costing also allowed for a comparison of various models of service delivery in different geographic regions from the perspective of service providers in Kenya. In addition, the analysis was also conducted to understand factors that indicate why some individuals place a greater value on PrEP than others, using a contingent valuation technique. Methods: Data collection was completed between November 2017 and September 2018. Costing data was collected from 44 Kenyan health facilities, consisting of 23 public facilities, 5 private facilities and 16 drop-in centers (DICEs) through a cross-sectional survey in ten counties. Financial and programmatic data were collected from financial and asset records and through interviewer administered questionnaires. The costs associated with PrEP provision were calculated using an ingredients-based costing approach which involved identification and costing of all the economic inputs (both direct and indirect) used in PrEP service delivery. In addition, a contingent valuation study was conducted at the same 44 facilities to understand factors that reveal why some individuals place a greater value on PrEP than others. Interviews were conducted with 2,258 individuals (1,940 current PrEP clients and 318 non-PrEP clients). A contingent valuation method using a "payment card approach" was used to determine the maximum willingness to pay (WTP) of respondents regarding obtaining access to oral PrEP services. Results: The weighted cost of providing PrEP was $253 per person year, ranging from $217 at health centers to $283 at dispensaries. Drop-in centers (DICEs), which served about two-thirds of the client volume at surveyed facilities, had a unit cost of $276. The unit cost was highest for facilities targeting MSM ($355), while it was lowest for those targeting FSW ($248). The unit cost for facilities targeting AGYW was $323 per person year. The largest percentage of costs were attributable to personnel (58.5%), followed by the cost of drugs, which represented 25% of all costs. The median WTP for PrEP was $2 per month (mean was $4.07 per month). This covers only one-third of the monthly cost of the medication (approximately $6 per month) and less than 10% of the full cost of delivering PrEP ($21 per month). A sizable proportion of current clients (27%) were unwilling to pay anything for PrEP. Certain populations put a higher value on PrEP services, including: FSW and MSM, Muslims, individuals with higher education, persons between the ages of 20 and 35, and households with a higher income and expenditures. Discussion: This is the most recent and comprehensive study on the cost of PrEP delivery in Kenya. These results will be used in determining resource requirements and for resource mobilization to facilitate sustainable PrEP scale-up in Kenya and beyond. This contingent valuation study does have important implications for Kenya's PrEP program. First, it indicates that some populations are more motivated to adopt oral PrEP, as indicated by their higher WTP for the service. MSM and FSW, for example, placed a higher value on PrEP than AGYW. Higher educated individuals, in turn, put a much higher value on PrEP than those with less education (which may also reflect the higher "ability to pay" among those with more education). This suggests that any attempt to increase demand or improve PrEP continuation should consider these differences in client populations. Cost recovery from existing PrEP clients would have potentially negative consequences for uptake and continuation.

Methods and Assumptions for Estimating Key HIV Indicators in the UNAIDS Annual Estimates Process.

BACKGROUND: Each year UNAIDS supports national teams to estimate key HIV indicators using their latest data. These estimates are produced using a collection of models and software tools. This paper describes the demographic and HIV projection models used in this process. METHODS: The demographic model (DemProj) projects the population by sex and single age for each year of the estimate. This information is fed into the HIV model (AIDS Impact Model) to estimate key HIV indicators. The model uses program, survey and surveillance data along with incidence trends produced through 1 of several separate models, to estimate new HIV infections, HIV-related deaths, and the population living with HIV by sex, age, CD4 category, and treatment status. RESULTS: These models allow the annual production of estimates of key HIV indicators including uncertainty intervals. This information is used to track progress toward national and global goals and to develop national strategic plans, Global Fund applications and PEPFAR country operational plans. CONCLUSIONS: Under the guidance of the UNAIDS Reference Group on Estimates, Modeling and Projections, these models are updated on a regular basis in response to evolving programmatic needs, new data, and analyses. This process of continuous review and improvement has led to mature models that make the best use of available data to provide estimates of indicators important to monitoring progress and developing future plans.

New HIV Infections Among Key Populations and Their Partners in 2010 and 2022, by World Region: A Multisources Estimation.

BACKGROUND: Previously, The Joint United Nations Programme on HIV/AIDS estimated proportions of adult new HIV infections among key populations (KPs) in the last calendar year, globally and in 8 regions. We refined and updated these, for 2010 and 2022, using country-level trend models informed by national data. METHODS: Infections among 15-49 year olds were estimated for sex workers (SWs), male clients of female SW, men who have sex with men (MSM), people who inject drugs (PWID), transgender women (TGW), and non-KP sex partners of these groups. Transmission models used were Goals (71 countries), AIDS Epidemic Model (13 Asian countries), Optima (9 European and Central Asian countries), and Thembisa (South Africa). Statistical Estimation and Projection Package fits were used for 15 countries. For 40 countries, new infections in 1 or more KPs were approximated from first-time diagnoses by the mode of transmission. Infection proportions among nonclient partners came from Goals, Optima, AIDS Epidemic Model, and Thembisa. For remaining countries and groups not represented in models, median proportions by KP were extrapolated from countries modeled within the same region. RESULTS: Across 172 countries, estimated proportions of new adult infections in 2010 and 2022 were both 7.7% for SW, 11% and 20% for MSM, 0.72% and 1.1% for TGW, 6.8% and 8.0% for PWID, 12% and 10% for clients, and 5.3% and 8.2% for nonclient partners. In sub-Saharan Africa, proportions of new HIV infections decreased among SW, clients, and non-KP partners but increased for PWID; elsewhere these groups' 2010-to-2022 differences were opposite. For MSM and TGW, the proportions increased across all regions. CONCLUSIONS: KPs continue to have disproportionately high HIV incidence.

Triangulation of Routine Antenatal HIV Prevalence Data and Adjusted HIV Estimates in Mozambique.

BACKGROUND: Routine health system data are central to monitoring HIV trends. In Mozambique, the reported number of women receiving antenatal care (ANC) and antiretroviral therapy for prevention of mother-to-child transmission (PMTCT) has exceeded the Spectrum-estimated number of pregnant women since 2017. In some provinces, reported HIV prevalence in pregnant women has declined faster than epidemiologically plausible. We hypothesized that these issues are linked and caused by programmatic overenumeration of HIV-negative pregnant women at ANC. METHODS: We triangulated program-reported ANC client numbers with survey-based fertility estimates and facility birth data adjusted for the proportion of facility births. We used survey-reported ANC attendance to produce adjusted time series of HIV prevalence in pregnant women, adjusted for hypothesized program double counting. We calibrated the Spectrum HIV estimation models to adjusted HIV prevalence data to produce adjusted adult and pediatric HIV estimates. RESULTS: ANC client numbers were not consistent with facility birth data or modeled population estimates indicating ANC data quality issues in all provinces. Adjusted provincial ANC HIV prevalence in 2021 was median 45% [interquartile range 35%-52% or 2.3 percentage points (interquartile range 2.5-3.5)] higher than reported HIV prevalence. In 2021, calibrating to adjusted antenatal HIV prevalence lowered PMTCT coverage to less than 100% in most provinces and increased the modeled number of new child infections by 35%. The adjusted results better reconciled adult and pediatric antiretroviral treatment coverage and antenatal HIV prevalence with regional fertility estimates. CONCLUSIONS: Adjusting HIV prevalence in pregnant women using nationally representative household survey data on ANC attendance produced estimates more consistent with surveillance data. The number of children living with HIV in Mozambique has been substantially underestimated because of biased routine ANC prevalence. Renewed focus on HIV surveillance among pregnant women would improve PMTCT coverage and pediatric HIV estimates.

Differences in Breastfeeding Duration by Maternal HIV Status: A Pooled Analysis of Nationally Representative Surveys in Sub-Saharan Africa.

BACKGROUND: Breastfeeding improves child survival but is a source of mother-to-child HIV transmission among women with unsuppressed HIV infection. Estimated HIV incidence in children is sensitive to breastfeeding duration among mothers living with HIV (MLHIV). Breastfeeding duration may vary according to maternal HIV status. SETTING: Sub-Saharan Africa. METHODS: We analyzed pooled data from nationally representative household surveys conducted during 2003-2019 that included HIV testing and elicited breastfeeding practices. We fitted survival models of breastfeeding duration by country, year, and maternal HIV status for 4 sub-Saharan African regions (Eastern, Central, Southern, and Western). RESULTS: Data were obtained from 65 surveys in 31 countries. In 2010, breastfeeding in the first month of life ("initial breastfeeding") among MLHIV ranged from 69.1% (95% credible interval: 68-79.9) in Southern Africa to 93.4% (92.7-98.0) in Western Africa. Median breastfeeding duration among MLHIV was the shortest in Southern Africa at 15.6 (14.2-16.3) months and the longest in Eastern Africa at 22.0 (21.7-22.5) months. By comparison, HIV-negative mothers were more likely to breastfeed initially (91.0%-98.7% across regions) and for longer duration (median 18.3-24.6 months across regions). Initial breastfeeding and median breastfeeding duration decreased during 2005-2015 in most regions and did not increase in any region regardless of maternal HIV status. CONCLUSIONS: MLHIV in sub-Saharan Africa are less likely to breastfeed initially and stop breastfeeding sooner than HIV-negative mothers. Since 2020, UNAIDS-supported HIV estimates have accounted for this shorter breastfeeding exposure among HIV-exposed children. MLHIV need support to enable optimal breastfeeding practices and to adhere to antiretroviral therapy for HIV treatment and prevention of postnatal mother-to-child transmission.

Estimation of Improvements in Mortality in Spectrum Among Adults With HIV Receiving Antiretroviral Therapy in High-Income Countries.

INTRODUCTION: Mortality rates for people living with HIV (PLHIV) on antiretroviral therapy (ART) in high-income countries continue to decline. We compared mortality rates among PLHIV on ART in Europe for 2016-2020 with Spectrum's estimates. METHODS: The AIDS Impact Module in Spectrum is a compartmental HIV epidemic model coupled with a demographic population projection model. We used national Spectrum projections developed for the 2022 HIV estimates round to calculate mortality rates among PLHIV on ART, adjusting to the age/country distribution of PLHIV starting ART from 1996 to 2020 in the Antiretroviral Therapy Cohort Collaboration (ART-CC)'s European cohorts. RESULTS: In the ART-CC, 11,504 of 162,835 PLHIV died. Between 1996-1999 and 2016-2020, AIDS-related mortality in the ART-CC decreased from 8.8 (95% CI: 7.6 to 10.1) to 1.0 (0.9-1.2) and from 5.9 (4.4-8.1) to 1.1 (0.9-1.4) deaths per 1000 person-years among men and women, respectively. Non-AIDS-related mortality decreased from 9.1 (7.9-10.5) to 6.1 (5.8-6.5) and from 7.0 (5.2-9.3) to 4.8 (4.3-5.2) deaths per 1000 person-years among men and women, respectively. Adjusted all-cause mortality rates in Spectrum among men were near ART-CC estimates for 2016-2020 (Spectrum: 7.02-7.47 deaths per 1000 person-years) but approximately 20% lower in women (Spectrum: 4.66-4.70). Adjusted excess mortality rates in Spectrum were 2.5-fold higher in women and 3.1-3.4-fold higher in men in comparison to the ART-CC's AIDS-specific mortality rates. DISCUSSION: Spectrum's all-cause mortality estimates among PLHIV are consistent with age/country-controlled mortality observed in ART-CC, with some underestimation of mortality among women. Comparing results suggest that 60%-70% of excess deaths among PLHIV on ART in Spectrum are from non-AIDS causes.

Modeling the epidemiological impact of the UNAIDS 2025 targets to end AIDS as a public health threat by 2030.

BACKGROUND: UNAIDS has established new program targets for 2025 to achieve the goal of eliminating AIDS as a public health threat by 2030. This study reports on efforts to use mathematical models to estimate the impact of achieving those targets. METHODS AND FINDINGS: We simulated the impact of achieving the targets at country level using the Goals model, a mathematical simulation model of HIV epidemic dynamics that includes the impact of prevention and treatment interventions. For 77 high-burden countries, we fit the model to surveillance and survey data for 1970 to 2020 and then projected the impact of achieving the targets for the period 2019 to 2030. Results from these 77 countries were extrapolated to produce estimates for 96 others. Goals model results were checked by comparing against projections done with the Optima HIV model and the AIDS Epidemic Model (AEM) for selected countries. We included estimates of the impact of societal enablers (access to justice and law reform, stigma and discrimination elimination, and gender equality) and the impact of Coronavirus Disease 2019 (COVID-19). Results show that achieving the 2025 targets would reduce new annual infections by 83% (71% to 86% across regions) and AIDS-related deaths by 78% (67% to 81% across regions) by 2025 compared to 2010. Lack of progress on societal enablers could endanger these achievements and result in as many as 2.6 million (44%) cumulative additional new HIV infections and 440,000 (54%) more AIDS-related deaths between 2020 and 2030 compared to full achievement of all targets. COVID-19-related disruptions could increase new HIV infections and AIDS-related deaths by 10% in the next 2 years, but targets could still be achieved by 2025. Study limitations include the reliance on self-reports for most data on behaviors, the use of intervention effect sizes from published studies that may overstate intervention impacts outside of controlled study settings, and the use of proxy countries to estimate the impact in countries with fewer than 4,000 annual HIV infections. CONCLUSIONS: The new targets for 2025 build on the progress made since 2010 and represent ambitious short-term goals. Achieving these targets would bring us close to the goals of reducing new HIV infections and AIDS-related deaths by 90% between 2010 and 2030. By 2025, global new infections and AIDS deaths would drop to 4.4 and 3.9 per 100,000 population, and the number of people living with HIV (PLHIV) would be declining. There would be 32 million people on treatment, and they would need continuing support for their lifetime. Incidence for the total global population would be below 0.15% everywhere. The number of PLHIV would start declining by 2023.

Updates to Spectrum's case surveillance and vital registration tool for HIV estimates and projections.

INTRODUCTION: The Case Surveillance and Vital Registration (CSAVR) model within Spectrum estimates HIV incidence trends from surveillance data on numbers of new HIV diagnoses and HIV-related deaths. This article describes developments of the CSAVR tool to more flexibly model diagnosis rates over time, estimate incidence patterns by sex and age group and by key population group. METHODS: We modelled HIV diagnosis rate trends as a mixture of three factors, including temporal and opportunistic infection components. The tool was expanded to estimate incidence rate ratios by sex and age for countries with disaggregated reporting of new HIV diagnoses and AIDS deaths, and to account for information on key populations such as men who have sex with men (MSM), males who inject drugs (MWID), female sex workers (FSW) and females who inject drugs (FWID). We used a Bayesian framework to calibrate the tool in 71 high-income or low-HIV burden countries. RESULTS: Across countries, an estimated median 89% (interquartile range [IQR]: 78%-96%) of HIV-positive adults knew their status in 2019. Mean CD4 counts at diagnosis were stable over time, with a median of 456 cells/µl (IQR: 391-508) across countries in 2019. In European countries reporting new HIV diagnoses among key populations, median estimated proportions of males that are MSM and MWID was 1.3% (IQR: 0.9%-2.0%) and 0.56% (IQR: 0.51%-0.64%), respectively. The median estimated proportions of females that are FSW and FWID were 0.36% (IQR: 0.27%-0.45%) and 0.14 (IQR: 0.13%-0.15%), respectively. HIV incidence per 100 person-years increased among MSM, with median estimates reaching 0.43 (IQR: 0.29-1.73) in 2019, but remained stable in MWID, FSW and FWID, at around 0.12 (IQR: 0.04-1.9), 0.09 (IQR: 0.06-0.69) and 0.13% (IQR: 0.08%-0.91%) in 2019, respectively. Knowledge of HIV status among HIV-positive adults gradually increased since the early 1990s to exceed 75% in more than 75% of countries in 2019 among each key population. CONCLUSIONS: CSAVR offers an approach to using routine surveillance and vital registration data to estimate and project trends in both HIV incidence and knowledge of HIV status.

Disease progression and mortality with untreated HIV infection: evidence synthesis of HIV seroconverter cohorts, antiretroviral treatment clinical cohorts and population-based survey data.

INTRODUCTION: Model-based estimates of key HIV indicators depend on past epidemic trends that are derived based on assumptions about HIV disease progression and mortality in the absence of antiretroviral treatment (ART). Population-based HIV Impact Assessment (PHIA) household surveys conducted between 2015 and 2018 found substantial numbers of respondents living with untreated HIV infection. CD4 cell counts measured in these individuals provide novel information to estimate HIV disease progression and mortality rates off ART. METHODS: We used Bayesian multi-parameter evidence synthesis to combine data on (1) cross-sectional CD4 cell counts among untreated adults living with HIV from 10 PHIA surveys, (2) survival after HIV seroconversion in East African seroconverter cohorts, (3) post-seroconversion CD4 counts and (4) mortality rates by CD4 count predominantly from European, North American and Australian seroconverter cohorts. We used incremental mixture importance sampling to estimate HIV natural history and ART uptake parameters used in the Spectrum software. We validated modelled trends in CD4 count at ART initiation against ART initiator cohorts in sub-Saharan Africa. RESULTS: Median untreated HIV survival decreased with increasing age at seroconversion, from 12.5 years [95% credible interval (CrI): 12.1-12.7] at ages 15-24 to 7.2 years (95% CrI: 7.1-7.7) at ages 45-54. Older age was associated with lower initial CD4 counts, faster CD4 count decline and higher HIV-related mortality rates. Our estimates suggested a weaker association between ART uptake and HIV-related mortality rates than previously assumed in Spectrum. Modelled CD4 counts in untreated people living with HIV matched recent household survey data well, though some intercountry variation in frequencies of CD4 counts above 500 cells/mm3 was not explained. Trends in CD4 counts at ART initiation were comparable to data from ART initiator cohorts. An alternate model that stratified progression and mortality rates by sex did not improve model fit appreciably. CONCLUSIONS: Synthesis of multiple data sources results in similar overall survival as previous Spectrum parameter assumptions but implies more rapid progression and longer survival in lower CD4 categories. New natural history parameter values improve consistency of model estimates with recent cross-sectional CD4 data and trends in CD4 counts at ART initiation.

Updates to the Spectrum/AIM model for the UNAIDS 2020 HIV estimates.

INTRODUCTION: The Spectrum/AIM model is used by national HIV programs and UNAIDS to prepare annual estimates of key HIV indicators. This article describes key updates to paediatric and adult models for the 2021 round of HIV estimates. METHODS: Potential updates to Spectrum arise due to newly available data, new analyses of existing data, and the need for new issues to be addressed. Updates are guided by experts through the UNAIDS Reference Group on Estimates, Modelling and Projections. Changes are tested and assessed for impact before being accepted into the final model. RESULTS: Spectrum tracks children living with HIV by CD4% for ages 0-4 and CD4 count for ages 5-14. Data from IeDEA treatment sites have been used to map the transition from CD4% to CD4 count at age 5. Breastfeeding patterns in sub-Saharan Africa have been updated with the latest survey data and estimates of continuation on antiretroviral therapy (ART) with breastfeeding have been revised based on recent studies. Model assumptions about the CD4 counts of people who drop out of ART have been revised to account for CD4 count increases while on treatment. If available, monthly data on numbers on ART can now be used to estimate the effects of COVID-19-related disruptions during 2020. CONCLUSIONS: These changes are intended to provide more accurate estimates of HIV burden. The effects of these changes on paediatric indicators are small except in countries with new surveys that might have updated patterns of breastfeeding. Changes to the adult model have little effect on total new infections. AIDS-related deaths will be somewhat lower in countries that have data on ART drop out but might be increased by HIV care disruptions due to COVID-19. The updated model uses newly available data to improve the estimation of paediatric and adult HIV indicators.

Naomi: a new modelling tool for estimating HIV epidemic indicators at the district level in sub-Saharan Africa.

INTRODUCTION: HIV planning requires granular estimates for the number of people living with HIV (PLHIV), antiretroviral treatment (ART) coverage and unmet need, and new HIV infections by district, or equivalent subnational administrative level. We developed a Bayesian small-area estimation model, called Naomi, to estimate these quantities stratified by subnational administrative units, sex, and five-year age groups. METHODS: Small-area regressions for HIV prevalence, ART coverage and HIV incidence were jointly calibrated using subnational household survey data on all three indicators, routine antenatal service delivery data on HIV prevalence and ART coverage among pregnant women, and service delivery data on the number of PLHIV receiving ART. Incidence was modelled by district-level HIV prevalence and ART coverage. Model outputs of counts and rates for each indicator were aggregated to multiple geographic and demographic stratifications of interest. The model was estimated in an empirical Bayes framework, furnishing probabilistic uncertainty ranges for all output indicators. Example results were presented using data from Malawi during 2016-2018. RESULTS: Adult HIV prevalence in September 2018 ranged from 3.2% to 17.1% across Malawi's districts and was higher in southern districts and in metropolitan areas. ART coverage was more homogenous, ranging from 75% to 82%. The largest number of PLHIV was among ages 35 to 39 for both women and men, while the most untreated PLHIV were among ages 25 to 29 for women and 30 to 34 for men. Relative uncertainty was larger for the untreated PLHIV than the number on ART or total PLHIV. Among clients receiving ART at facilities in Lilongwe city, an estimated 71% (95% CI, 61% to 79%) resided in Lilongwe city, 20% (14% to 27%) in Lilongwe district outside the metropolis, and 9% (6% to 12%) in neighbouring Dowa district. Thirty-eight percent (26% to 50%) of Lilongwe rural residents and 39% (27% to 50%) of Dowa residents received treatment at facilities in Lilongwe city. CONCLUSIONS: The Naomi model synthesizes multiple subnational data sources to furnish estimates of key indicators for HIV programme planning, resource allocation, and target setting. Further model development to meet evolving HIV policy priorities and programme need should be accompanied by continued strengthening and understanding of routine health system data.

Strategic options for syphilis control in Papua New Guinea- impact and cost-effectiveness projections using the syphilis interventions towards elimination (SITE) model.

OBJECTIVES: Papua New Guinea (PNG) has among the highest rates of sexually transmitted infections (STIs) globally and is committed to reducing their incidence. The Syphilis Interventions Towards Elimination (SITE) model was used to explore the expected impact and cost of alternative syphilis intervention scale-up scenarios. METHODS: SITE is a dynamical model of syphilis transmission among adults 15-49 years. Individuals are divided into nine groups based on sexual behaviour and into six stages of infection. The model was calibrated to PNG using data from routine surveillance, bio-behavioural surveys, research studies and program records. Inputs included syphilis prevalence, risk behaviours, intervention coverage and service delivery unit costs. Scenarios compared different interventions (clinical treatment, contact tracing, syphilis screening, and condom promotion) for incidence and cost per infection averted over 2021-2030. RESULTS: Increasing treatment coverage of symptomatic primary/secondary-stage syphilis cases from 25-35% in 2020 to 60% from 2023 onwards reduced estimated incidence over 2021-2030 by 55%, compared to a scenario assuming constant coverage at 2019-2020 levels. The introduction of contact tracing in 2020, assuming 0.4 contacts per symptomatic person treated, reduced incidence over 2021-2030 by 10%. Increasing screening coverage by 20-30 percentage points from the 2019-2020 level reduced incidence over 2021-2030 by 3-16% depending on the target population. Scaling-up clinical, symptom-driven treatment and contact tracing had the lowest cost per infection averted, followed by condom promotion and periodic screening of female sex workers and men who have sex with men. CONCLUSIONS: PNG could considerably reduce its syphilis burden by scaling-up clinical treatment and contact tracing alongside targeted behavioural risk reduction interventions. SITE is a useful tool countries can apply to inform national STI programming and resource allocation.

Potential effects of disruption to HIV programmes in sub-Saharan Africa caused by COVID-19: results from multiple mathematical models.

BACKGROUND: The COVID-19 pandemic could lead to disruptions to provision of HIV services for people living with HIV and those at risk of acquiring HIV in sub-Saharan Africa, where UNAIDS estimated that more than two-thirds of the approximately 38 million people living with HIV resided in 2018. We aimed to predict the potential effects of such disruptions on HIV-related deaths and new infections in sub-Saharan Africa. METHODS: In this modelling study, we used five well described models of HIV epidemics (Goals, Optima HIV, HIV Synthesis, an Imperial College London model, and Epidemiological MODeling software [EMOD]) to estimate the effect of various potential disruptions to HIV prevention, testing, and treatment services on HIV-related deaths and new infections in sub-Saharan Africa lasting 6 months over 1 year from April 1, 2020. We considered scenarios in which disruptions affected 20%, 50%, and 100% of the population. FINDINGS: A 6-month interruption of supply of antiretroviral therapy (ART) drugs across 50% of the population of people living with HIV who are on treatment would be expected to lead to a 1·63 times (median across models; range 1·39-1·87) increase in HIV-related deaths over a 1-year period compared with no disruption. In sub-Saharan Africa, this increase amounts to a median excess of HIV deaths, across all model estimates, of 296 000 (range 229 023-420 000) if such a high level of disruption occurred. Interruption of ART would increase mother-to-child transmission of HIV by approximately 1·6 times. Although an interruption in the supply of ART drugs would have the largest impact of any potential disruptions, effects of poorer clinical care due to overstretched health facilities, interruptions of supply of other drugs such as co-trimoxazole, and suspension of HIV testing would all have a substantial effect on population-level mortality (up to a 1·06 times increase in HIV-related deaths over a 1-year period due to disruptions affecting 50% of the population compared with no disruption). Interruption to condom supplies and peer education would make populations more susceptible to increases in HIV incidence, although physical distancing measures could lead to reductions in risky sexual behaviour (up to 1·19 times increase in new HIV infections over a 1-year period if 50% of people are affected). INTERPRETATION: During the COVID-19 pandemic, the primary priority for governments, donors, suppliers, and communities should focus on maintaining uninterrupted supply of ART drugs for people with HIV to avoid additional HIV-related deaths. The provision of other HIV prevention measures is also important to prevent any increase in HIV incidence. FUNDING: Bill & Melinda Gates Foundation.

Modelling impact and cost-effectiveness of oral pre-exposure prophylaxis in 13 low-resource countries.

INTRODUCTION: Oral pre-exposure prophylaxis (PrEP) provision is a priority intervention for high HIV prevalence settings and populations at substantial risk of HIV acquisition. This mathematical modelling analysis estimated the impact, cost and cost-effectiveness of scaling up oral PrEP in 13 countries. METHODS: We projected the impact and cost-effectiveness of oral PrEP between 2018 and 2030 using a combination of the Incidence Patterns Model and the Goals model. We created four PrEP rollout scenarios involving three priority populations-female sex workers (FSWs), serodiscordant couples (SDCs) and adolescent girls and young women (AGYW)-both with and without geographic prioritization. We applied the model to 13 countries (Eswatini, Ethiopia, Haiti, Kenya, Lesotho, Mozambique, Namibia, Nigeria, Tanzania, Uganda, Zambia and Zimbabwe). The base case assumed achievement of the Joint United Nations Programme on HIV/AIDS 90-90-90 antiretroviral therapy targets, 90% male circumcision coverage by 2020 and 90% efficacy and adherence levels for oral PrEP. RESULTS: In the scenarios we examined, oral PrEP averted 3% to 8% of HIV infections across the 13 countries between 2018 and 2030. For all but three countries, more than 50% of the HIV infections averted by oral PrEP in the scenarios we examined could be obtained by rollout to FSWs and SDCs alone. For several countries, expanding oral PrEP to include medium-risk AGYW in all regions greatly increased the impact. The efficiency and impact benefits of geographic prioritization of rollout to AGYW varied across countries. Variations in cost-effectiveness across countries reflected differences in HIV incidence and expected variations in unit cost. For most countries, rolling out oral PrEP to FSWs, SDCs and geographically prioritized AGYW was not projected to have a substantial impact on the supply chain for antiretroviral drugs. CONCLUSIONS: These modelling results can inform prioritization, target-setting and other decisions related to oral PrEP scale-up within combination prevention programmes. We caution against extensive use given limitations in cost data and implementation approaches. This analysis highlights some of the immediate challenges facing countries-for example, trade-offs between overall impact and cost-effectiveness-and emphasizes the need to improve data availability and risk assessment tools to help countries make informed decisions.

The Estimation and Projection Package Age-Sex Model and the r-hybrid model: new tools for estimating HIV incidence trends in sub-Saharan Africa.

OBJECTIVES: Improve models for estimating HIV epidemic trends in sub-Saharan Africa (SSA). DESIGN: Mathematical epidemic model fit to national HIV survey and ANC sentinel surveillance (ANC-SS) data. METHODS: We modified EPP to incorporate age and sex stratification (EPP-ASM) to more accurately capture the shifting demographics of maturing HIV epidemics. Secondly, we developed a new functional form for the HIV transmission rate, termed 'r-hybrid', which combines a four-parameter logistic function for the initial epidemic growth, peak, and decline followed by a first-order random walk for recent trends after epidemic stabilization. We fitted the r-hybrid model along with previously developed r-spline and r-trend models to HIV prevalence data from household surveys and ANC-SS in 177 regions in 34 SSA countries. We used leave-one-out cross validation with household survey HIV prevalence to compare model predictions. RESULTS: The r-hybrid and r-spline models typically provided similar HIV prevalence trends, but sometimes qualitatively different assessments of recent incidence trends because of different structural assumptions about the HIV transmission rate. The r-hybrid model had the lowest average continuous ranked probability score, indicating the best model predictions. Coverage of 95% posterior predictive intervals was 91.5% for the r-hybrid model, versus 87.2 and 85.5% for r-spline and r-trend, respectively. CONCLUSION: The EPP-ASM and r-hybrid models improve consistency of EPP and Spectrum, improve the epidemiological assumptions underpinning recent HIV incidence estimates, and improve estimates and short-term projections of HIV prevalence trends. Countries that use general population survey and ANC-SS data to estimate HIV epidemic trends should consider using these tools.

Updates to the Spectrum/AIM model for estimating key HIV indicators at national and subnational levels.

BACKGROUND: The Spectrum/AIM model is used by national programs and UNAIDS to prepare annual estimates of the status of the HIV epidemic in 170 countries. The model and assumptions are updated regularly under the guidance of the UNAIDS Reference Group on Estimates, Modelling and Projections in response to new data, studies and program needs. This article describes the most recent updates for the 2018 round of estimates. METHODS: New data on AIDS-related mortality from Europe and Brazil have been used to update mortality rates of those not on antiretroviral therapy (ART). Household survey data and new studies of pregnant women, mothers, and children have been used to improve estimates of the number of HIV-positive pregnant and breastfeeding women and pediatric ART initiation. New tools to estimate geographic variation in HIV prevalence have been used to prepare district estimates of key indicators. RESULTS: The 2018 version of Spectrum includes: new estimates of non-ART AIDS-related mortality by CD4 count that depend on ART coverage; a procedure to estimate country-specific patterns of HIV incidence by age by fitting to prevalence by age from household surveys; an updated estimate of postpartum transmission with ART started before pregnancy of 0.023% per month; an updated estimate of retention on treatment at delivery of 80% for all women on ART; a somewhat older pattern of ART initiation by age that has 26% of new pediatric patients initiating ART at 10-14 years of age, 18% at 2-4 years of age, and 26% at 5-9 years of age; and a new tool for estimating key HIV indicators at the district level. CONCLUSION: The new methods and data implemented in the 2018 version of Spectrum allow national programs more flexibility in describing their programs and are intended to improve the estimates of adult mortality and pediatric HIV indicators.

Estimating and projecting the number of new HIV diagnoses and incidence in Spectrum's case surveillance and vital registration tool.

OBJECTIVE: The Joint United Nations Programme on HIV/AIDS-supported Spectrum software package is used by most countries worldwide to monitor the HIV epidemic. In Spectrum, HIV incidence trends among adults (aged 15-49 years) are derived by either fitting to seroprevalence surveillance and survey data or generating curves consistent with case surveillance and vital registration data, such as historical trends in the number of newly diagnosed infections or AIDS-related deaths. This article describes development and application of the case surveillance and vital registration (CSAVR) tool for the 2019 estimate round. METHODS: Incidence in CSAVR is either estimated directly using single logistic, double logistic, or spline functions, or indirectly via the 'r-logistic' model, which represents the (log-transformed) per-capita transmission rate using a logistic function. The propensity to get diagnosed is assumed to be monotonic, following a Gamma cumulative distribution function and proportional to mortality as a function of time since infection. Model parameters are estimated from a combination of historical surveillance data on newly reported HIV cases, mean CD4 at HIV diagnosis and estimates of AIDS-related deaths from vital registration systems. Bayesian calibration is used to identify the best fitting incidence trend and uncertainty bounds. RESULTS: We used CSAVR to estimate HIV incidence, number of new diagnoses, mean CD4 at diagnosis and the proportion undiagnosed in 31 European, Latin American, Middle Eastern, and Asian-Pacific countries. The spline model appeared to provide the best fit in most countries (45%), followed by the r-logistic (25%), double logistic (25%), and single logistic models. The proportion of HIV-positive people who knew their status increased from about 0.31 [interquartile range (IQR): 0.10-0.45] in 1990 to about 0.77 (IQR: 0.50-0.89) in 2017. The mean CD4 at diagnosis appeared to be stable, at around 410 cells/µl (IQR: 224-567) in 1990 and 373 cells/µl (IQR: 174-475) by 2017. CONCLUSION: Robust case surveillance and vital registration data are routinely available in many middle-income and high-income countries while HIV seroprevalence surveillance and survey data may be scarce. In these countries, CSAVR offers a simpler, improved approach to estimating and projecting trends in both HIV incidence and knowledge of HIV status.

Drug resistance from preferred antiretroviral regimens for HIV infection in South Africa: A modeling study.

BACKGROUND: Tenofovir-containing regimens comprise the preferred first-line antiretroviral therapy (ART) in many countries including South Africa, where utilization of second-line regimens is limited. Considerable HIV drug resistance has occurred among persons failing tenofovir-containing first-line ART. We evaluated drug resistance at the population level using mathematical modeling. SETTING: Heterosexual HIV epidemic in KwaZulu-Natal, South Africa. METHODS: We constructed a stochastic individual-based model and simulated scenarios of ART implementation, either CD4-based (threshold < 500 cells/mL) or Fast-track (81% coverage by 2020), with consideration of major drug-associated mutations (M184V, K65R and non-nucleoside reverse transcriptase inhibitor (NNRTI)). Using base case and uncertainty analyses, we assessed (majority) drug resistance levels. RESULTS: By 2030, the median total resistance (proportion of HIV-infected persons with drug resistance) is predicted to reach 31.4% (interquartile range (IQR): 16.5%-50.2%) with CD4-based ART, decreasing to 14.5% (IQR: 7.7%-25.8%) with Fast-track implementation. In both scenarios, we find comparably high prevalence (~80%) of acquired NNRTI-associated, M184V and K65R mutations. Over 48% of individuals with acquired resistance harbor dual, 44% triple and 7% just single drug mutations. Drug-resistant HIV is predicted to comprise 40% (IQR: 27%-50%) of incident infections, while 70% of prevalent transmitted resistance is NNRTI-associated. At 2018, the projected total resistance is 15% (IQR: 7.5%-25%), with 18% (IQR: 13%-24%) of incident infections from transmitted drug-resistant HIV. CONCLUSIONS: WHO-recommended preferred first-line ART could lead to substantial drug resistance. Effective surveillance of HIV drug resistance and utilization of second-line as well as alternative first-line regimens is crucial.

Evaluating the potential impact and cost-effectiveness of dapivirine vaginal ring pre-exposure prophylaxis for HIV prevention.

BACKGROUND: Expanded HIV prevention options are needed to increase uptake of HIV prevention among women, especially in generalized epidemics. As the dapivirine vaginal ring moves forward through regulatory review and open-label extension studies, the potential public health impact and cost-effectiveness of this new prevention method are not fully known. We used mathematical modeling to explore the impact and cost-effectiveness of the ring in different implementation scenarios alongside scale-up of other HIV prevention interventions. Given the knowledge gaps about key factors influencing the ring's implementation, including potential uptake and delivery costs, we engaged in a stakeholder consultation process to elicit plausible parameter ranges and explored scenarios to identify the possible range of impact, cost, and cost-effectiveness. METHODS AND FINDINGS: We used the Goals model to simulate scenarios of oral and ring pre-exposure prophylaxis (PrEP) implementation among female sex workers and among other women ≤21 years or >21 years with multiple male partners, in Kenya, South Africa, Uganda, and Zimbabwe. In these scenarios, we varied antiretroviral therapy (ART) coverage, dapivirine ring coverage and ring effectiveness (encompassing efficacy and adherence) by risk group. Following discussions with stakeholders, the maximum level of PrEP coverage (oral and/or ring) considered in each country was equal to modern contraception use minus condom use in the two age groups. We assessed results for 18 years, from 2018 to 2035. In South Africa, for example, the HIV infections averted by PrEP (ring plus oral PrEP) ranged from 310,000 under the highest-impact scenario (including ART held constant at 2017 levels, high ring coverage, and 85% ring effectiveness) to 55,000 under the lowest-impact scenario (including ART reaching the UNAIDS 90-90-90 targets by 2020, low ring coverage, and 30% ring effectiveness). This represented a range of 6.4% to 2.2% of new HIV infections averted. Given our assumptions, the addition of the ring results in 11% to 132% more impact than oral PrEP alone. The cost per HIV infection averted for the ring ranged from US$13,000 to US$121,000. CONCLUSIONS: This analysis offers a wide range of scenarios given the considerable uncertainty over ring uptake, consistency of use, and effectiveness, as well as HIV testing, prevention, and treatment use over the next two decades. This could help inform donors and implementers as they decide where to allocate resources in order to maximize the impact of the dapivirine ring in light of funding and implementation constraints. Better understanding of the cost and potential uptake of the intervention would improve our ability to estimate its cost-effectiveness and assess where it can have the most impact.

Dapivirine vaginal ring for HIV prevention: modelling health outcomes, drug resistance and cost-effectiveness.

INTRODUCTION: A vaginal ring containing dapivirine is effective for HIV prevention as pre-exposure prophylaxis (PrEP). We evaluated the potential epidemiological impact and cost-effectiveness of dapivirine vaginal ring PrEP among 22- to 45-year-old women in KwaZulu-Natal, South Africa. METHODS: Using mathematical modelling, we studied dapivirine vaginal ring PrEP implementation, either unprioritized, or prioritized based on HIV incidence (≥3% per year), age (22 to 29 years) or female sex worker status, alongside the implementation of voluntary medical male circumcision and antiretroviral therapy scaled-up to UNAIDS Fast-Track targets. Outcomes over the intervention (2019 to 2030) and lifetime horizons included cumulative HIV infections, life-years lived, costs and cost-effectiveness. We assessed the incremental cost-effectiveness ratios against the revealed willingness to pay ($500) and the standard (2017 per capita gross domestic product; $6161) cost-effectiveness thresholds for South Africa. RESULTS: Compared to a reference scenario without PrEP, implementation of dapivirine vaginal ring PrEP, assuming 56% effectiveness and covering 50% of 22 to 29-year-old or high-incidence women, prevented 10% or 11% of infections by 2030 respectively. Equivalent, unprioritized coverage (30%) prevented fewer infections (7%), whereas 50% coverage of female sex workers had the least impact (4%). Drug resistance attributable to PrEP was modest (2% to 4% of people living with drug-resistant HIV). Over the lifetime horizon, dapivirine PrEP implementation among female sex workers was cost-saving, whereas incidence-based PrEP cost $1898 per life-year gained, relative to PrEP among female sex workers and $989 versus the reference scenario. In a scenario of 37% PrEP effectiveness, PrEP had less impact, but prioritization to female sex workers remained cost-saving. In uncertainty analysis, female sex worker PrEP was consistently cost-saving; and over the lifetime horizon, PrEP cost less than $6161 per life-year gained in over 99% of simulations, whereas incidence- and age-based PrEP cost below $500 per life-year gained in 61% and 49% of simulations respectively. PrEP adherence and efficacy, and the effectiveness of antiretroviral therapy for HIV prevention, were the principal drivers of uncertainty in the cost-effectiveness of PrEP. CONCLUSIONS: Dapivirine vaginal ring PrEP would be cost-saving in KwaZulu-Natal if prioritized to female sex workers. PrEP's impact on HIV prevention would be increased, with potential affordability, if prioritized to women by age or incidence.